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Abstract 

We use the 1/M expansion for the CP^^^^ model to study the long-distance 
behaviour of the staggered spin susceptibility in the commensurate, two- 
dimensional quantum antiferromagnet at finite temperature. At M = oo this 
model possesses deconfined spin-1/2 bosonic spinons (Schwinger bosons), and 
the susceptibility has a branch cut along the imaginary k axis. We show that 
in all three scaling regimes at finite T, the interaction between spinons and 
gauge field fluctuations leads to divergent 1/M corrections near the branch 
cut. We identify the most divergent corrections to the susceptibility at each 
order in 1/M and explicitly show that the full static staggered susceptibility 
has a number of simple poles rather than a branch cut. We compare our 
results with the 1/A^ expansion for the 0{N) sigma-model. 
PACS: 67.50-b, 67.70-hn, 67.50Dg 
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I. INTRODUCTION 



Numerous studies of the two-dimensional (2D) quantum Heisenberg antiferromagnet 
(QHAF) undertaken in the last few years have significantly improved our understanding 
of the behavior of these systems at low temperatures There is no ordered state at any 

finite T (and hence, no phase transitions), but there are nevertheless three distinct low-T 
regimes depending on the relative values of the temperature and the coupling constant, g. 
These three regimes are [^: (i) renormalized-classical (RC) regime, where T is smaller than 
the spin stiffness in the ordered ground state, Qs] (n) quantum-disordered (QD) regime, 
when g is larger than the critical coupling for the T = disordering transition, and T is 
smaller than the gap. A, between the singlet ground state and the lowest excited state with 
S = 1; (Hi) quantum-critical (QC) regime, which lies between the other two, and in which 
temperature is the largest infrared cutoff, ksT ^ f?s, A. 

It is very likely that the low-energy physics of a 2D QHAF, at least in the vicinity of 
the disordering transition, is adequately described by the 0(3) nonlinear sigma- model (see 
Sec p. Some information about the properties of this model can be obtained from the 
Bethe-ansatz solution @J^, but most of the thermodynamic properties have been studied 
using several available perturbative techniques. The first perturbative approach to the 0(3) 
sigma-model was initiated many years ago by Polyakov In this approach, one departs 
from the ordered state at T = , and applies renormalization group theory which accounts 
for the effects of classical fluctuations at small but finite T. The expansion parameter 
for the RG studies is T/2TTgs- It has to be small, i.e., the system should be in the RC 
regime. The RG approach allows one to obtain a scale at which the fluctuation corrections 
to the spin-wave coupling constant become comparable to its bare value. This scale is then 
identified with the correlation length in a system. Calculations along these lines yielded 
e = A{hc/27rgs) exp{2Tr Qs/kBT) and S(g) = A^^ {kBT/2ngsfefiqO i- Here A is a 
constant {A ^ 0.34 0), S{q) is the static structure factor, A'o is the staggered magnetization, 
and f{x) tends to a finite value at a; = 0. 
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Another widely used approach is the expansion for the 0{N) sigma- model. The 
physical results are obtained in this approach by extending the perturbation series in 
to = 3. The advantage of the expansion is that it works in all three scaling regimes, 
and the point of departure is a disordered state at any finite temperature. The weak point of 
the theory is the absence of the physically motivated small parameter for the physical case 

= 3. In the RC regime, however, the series of 1/A^ terms can be explicitly summed up, and 
for the physical case of A^ = 3 one obtains exactly the same results as in the RG approach. At 
arbitrary A^, one finds ^ = AN{nc/kBT){{N - 2)kBT exp(27r^,/((A^ - 2)kBT) 
and S{q) = (fcBT/27rf?,)(^-i)/(^-2) ^'/^(gO (^3 = A, h{x) = f{x)). The results of 
1/N expansion for all three scaling regimes are collected in . 

Finally, the third approach is based on the spinon representation for spin operators 



(Schwinger boson theory). Read and Sachdev have shown explicitly [|T5[ that the low-energy 
limit of the Schwinger boson theory is described by a CP^ sigma-model for a two-component 
complex unit field (see below). To obtain this model, one has to express the unit vector field 
of the 0(3) sigma-model as a bilinear combination of two Bose fields {CP^ representation), 

"-a = 4</3^/3, (1-1) 

where are SU{2) indexes, a = x,y,z, and a" are the Pauli matrices. This represen- 
tation also introduces a U{1) gauge degree of freedom, because n remains invariant under 
the transformation of the bosonic fields, Za{'r,t) 2;„(r, t)e*'^*- Each z— field quantum 
carries S = 1/2 and is therefore a bosonic spinon. The condition = 1, however, imposes a 
local constraint zl,Za = 1, which implies that spinous can appear only in pairs. A mean-field 
version of the Schwinger-boson theory has recently received a lot of attention P). In the 
mean-field approximation, one reduces the on-site constraint to a constraint imposed on 
the averaged quantities, and decouples the term that is quartic in z in the spin Hamilto- 
nian. From the solution of the self-consistent equations in the RC regime one then obtains 
the spin correlation length, ^ ~ (ksT/ g^) exp{2ngs/kBT), and the static structure factor, 
S{q) ~ {kBT / 27!- QsY^^ f{q^), where f{x) has the same asymptotic behavior as f{x) [§. 



The temperature dependence of S{q) and the exponent in the expression for the correla- 
tion length are the same as in other approaches but the correlation length aquires an extra 
power of T in the prefactor. Because of this incorrect prefactor in ^, the validity of the 
Schwinger-boson mean-field theory has been questioned [^. More recently, however, the 
Schwinger-boson approach to a 2D antiferromagnet has been applied in a systematic way. 



by means of a controllable 1/M expansion |]TO|,|TI|. To generate this expansion, each z 



field was assumed to have M components rather than two (this changes the symmetry of 
the underlying sigma-model to CP^^~^). The 1/M computations have been performed for 
field-theory [|l^ and condensed-matter ||lT| applications. It has been shown that the extra 
power of T in the prefactor is indeed an artifact of the mean- field, M = oo, approximation: 
the correct prefactor is yi-^/M^ -^^ reduces to a constant for the physical case of M = 2. 
At the same time, the 1/M corrections to the static structure factor do not change the 
power of T, and the mean-field result for S{q) is therefore valid for all M. Clearly then, for 
physical spins, the n— field and the Schwinger-boson approaches yield the same temperature 
dependence of the observables in the RC regime, as they indeed should. One can therefore 
safely use any of these perturbative techniques. Notice however, that the series of regular 
1/M corrections are poorly convergent while the regular 1/A^ corrections are usually small. 
This makes the Schwinger boson approach less reliable for practical purposes than the 1/A^ 
expansion for the 0{N) sigma-model. 

There is, however, another discrepancy between the Schwinger-boson and the n— field 
approaches, which in our opinion has not been fully clarified in the literature. Namely, in the 
n— field approach, the staggered static spin susceptibility Xs{k,uj = 0) is proportional to the 
static n— field propagator which has a simple pole along the imaginary k axis, at A; = 
The residue of the pole is finite at = oo, and remains finite in the physical case of = 3. 
This result is valid in all three scaling regimes. On the other hand, in the Schwinger-boson 
formalism, the staggered spin susceptibility is a convolution of two spinon propagators. 
At the mean-field level, spinous behave as free particles, and elementary calculations show 
that the staggered susceptibility has only a branch cut singularity at A; = 2imo, where 



mo is the mass (inverse correlation length) of a Schwinger boson. Since the behavior near 
the singularity in 0) determines the long-distance properties of the spin correlators, the 
difference in the type of the singularity in the two models leads to different predictions about 
the long-distance behavior of the correlation function. Obviously, one of these predictions 
must be wrong. 

In this paper, we show that the branch cut behavior of 0) is also an artifact of 
the mean-field Schwinger-boson formalism. We will see that 1/M corrections are divergent 
near the branch cut, and eventually transform the branch cut into a simple pole. This 
phenomenon is closely related to the confinement of spinous in the CP^^^^ model in 2 + 



dimensions, first studied by Witten [|T2|. He found that massless gauge fluctuations in the 
(jpM-i ]2iodel give rise to a linear confining potential between spinous, and this yields a 
bound state with a mass, m, which is a nonanalytical function of 1/M: m = 2mo(l + 
0(1/M2/3)). This result was reproduced in a number of more recent papers (for a review 



see, e.g., Ref. ||T3[ and references therein). However, to the best of our knowledge, the effect 
of the 1/M corrections on the staggered susceptibility has not been studied in detail. The 
results of such study will be reported in this paper. Besides the RC regime, we will also 
study staggered susceptibility in the QD and QC regimes. 

Before we proceed to the description of our calculations, it is useful to specify which 1/M 
corrections are essential to our analysis. The point is that in the RC regime, there exists 
a self-energy correction to the spinon propagator of the form l/Mloglog(/c5T/mo). Since 
mo is itself exponential in T, the double logarithm in fact reduces to l/M\og{gs/kBT). A 
series of these logarithmic terms give rise to the above mentioned change in the power of 
temperature in the preexponential factor in ^ from T to J^i^^/m ggjQ^ assume that 

these corrections are already included into the expression for the Schwinger boson mass. We 
will therefore consider only regular 1/M corrections which, as we will show, are responsible 
for the confinement. 

The paper is organized as follows. In Sec. |T|, we will briefly review the large M expansion 
for the CP^^~^ sigma model and present M = oo results for the correlation length, spin 



susceptibility, and gauge field propagator. In Sec. |T|, we consider the static staggered 
susceptibility in the RC regime. We first compute the lowest-order 1/M corrections, select 
the most divergent ones, and then sum up the ladder series of divergent 1/M terms by 
reducing the problem to a Schrodinger equation. Discrete solutions of this Schrodinger 
equation will correspond to the poles in the staggered susceptibility. In Sec |IV| and |V| we 
report analogous calculations for QD and QC phases, respectively. Finally, in Sec ^ we 
state our conclusions and discuss open questions. 



II. THE SIGMA-MODEL 



Our starting point is the partition function for the 0(3) nonlinear sigma model in Eu- 
clidian space 



Z = J Dni5{n] - l)exp | j^J'^ dr j rf' 



(2.1) 



where i = x,y] I = 1,2, 3, and ^° and cq are the bare spin stiffness and spin-wave velocity. 
For simplicity, throughout the paper we choose the units where h = 1 and cq = 1. Vector n 
describes local staggered magnetization. In the CP^ representation (|1 . 1| ) , Z transforms into 

Z = f DzDz6{\z\^ - 1) exp\-2dl j^'^ dr f d^r \\d^z\^ - \zd^z 



(2.2) 

Here ^ = T,x,y. Introducing the Hubbard- Stratonovich vector field to decouple the 
quartic term, we obtain 

Z = j DzDzDA^5{\z\^ - I) exp|-2^° j^'^ '^'^ j ^^r|(9^ - iA^)^p| . (2.3) 

Now we generalize the doublet z to the M-component complex vector, rescale the z field to 
z —> \/Mz, and introduce the coupling constant g = M/2g^. Introducing then the constraint 
into the action in a standard way, we obtain for the partition function 



Z = J DzDzDA^DXexpl-^ J^^^ dr I dh |(9^ - iA^)^!^ + iA(|zp - M) 



(2.4) 
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This is the partition function for the CP^ ^ model [|14|. The action in ( p.4| ) is quadratic in 



the 2;, z fields, and we can therefore integrate them out. This generates an effective action for 
the and A fields, which contains M only as a prefactor 0,0]. At large M, this effective 



action can be well approximated by the quadratic fluctuations of and A around their 
saddle-point values [^] 

i < A >= < >= 0, (2.5) 

The value of the spinon mass, m, can be obtained in the 1/M expansion by solving the 
constraint equation to any given order in 1/M. The nonzero solution for m implies that 
the rotational symmetry is unbroken. A more detailed discussion of the derivation of the 



CP model can be found in ref. [15| 



At M = cxD, the theory is particularly simple, because it describes free massive spinous. 
The spinon Green's function Go(pi'^n) is given by 

GM-n) = -^^^T^^. (2.6) 
and the constraint equation reads 

^E. / 7^,G,{p,u;r.) = -. (2.7) 
Using the Pauli-Villars regularization of ultraviolet divergencies, we obtain 

i?-^-(-^-)-e-5). (-) 

where (yf^ = and A is the ultraviolet regulator. Depending on the values of g/gc 
and the temperature, the solutions of (p.8|) are (i) RC regime {qs > kBT): 

^0 = ksT exp —{AnQs/MT), where = (^/2) (1/(7 — l/s'c) is the renormalized stiffness; 
(a) QD regime (mo > ksT): ttiq = A + 0(exp — A/T), where A = 8TT{l/gc — l/g)', (Hi) QC 
regime (T > ^„ A): mo = efc^T, where 6 = 2 log[(V5 + l)/2]. 

A. Staggered susceptibility 

The static staggered susceptibility is defined in continuum limit as 
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Xik,0)5at = Y I " < 5,(r,r)5,(0,0) > expH(^ + Q) ■ r], (2.9) 

where is a small momentum, Q = (vr/a, vr/a), and A^s is the number of spins in the 
system. As each spin is a byproduct of 2;— fields, the physical susceptibility is related to the 
polarization operator of z, n(/c, 0) = Il{k) = TJ^lj j{d?q/^'n'^)G{q,uj)G{k + q,uj). In the 
RC and QC regimes, the relation between x{k, 0) and n(A;) can be obtained in the same 
way as in 0. We find 

x{k,0) = -^^U{k), (2.10) 

where A^o and Qs are the fully renormalized on-site magnetization and spin-stiffness at T = 0. 
Notice that there is a factor of 2 difference with the analogous expression for frustrated 
systems 0. In the QD regime, the rescaling factor between susceptibility and polarization 
operator can be related to the overall factor in the local susceptibility 0. However, in this 
regime, we will only obtain Il{k) up to an overall factor, so there is no need to discuss the 
exact relation between x{k) and Il{k) in the QD regime. 



B. Polarization operator 



We now present the expressions for the polarization operator at M = cxo. In the RC 
regime, the summation over frequency reduces to the u = term, and one obtains 

£p 1 1 



Hoik) = ksT J 



4vr2 (p-A?/2)2 + mg {p + k/2y + ml' 
The momentum integration yields P,|3|JT0[] 

ksT 



(2.11) 



no(fc) 



47r5.M2 - ml 



log 



52 — mn + 5 



mo 



(2.12) 



where we introduced 6^ = A;^/4 + mg. Clearly, Ilo{k) has a branch cut singularity at 5 = 0, 



i.e., at k"^ = —Ami. Near the singularity, we obtain 



no(A;) 



ksT 
8mo5 
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(2.13) 



The behavior of n(fc) near the branch cut determines the long-distance behavior of the 
spin susceptibihty in real space. Evaluating the Fourier transform of ( |2.13| ), we obtain 
x(r) ~ exp{—r/^)/r, where ^ = l/2mo at the mean-field level. At the same time, the 
susceptibility in the 0{N) sigma- model has a simple pole at fc^ = — and its long-distance 
behavior is ~ exp(— r/^)/-,/?. 

For the QD regime, the leading term in n(A;) is the T = piece. The frequency summa- 
tion is then replaced by the integration, and, to leading order in S, one obtains 

We see that the polarization operator has only a weak logarithmical singularity. In real 
space, the singularity in Ilo{k) gives rise to xi'f') ~ exp(— 2mor)/r^, which is again different 
from the mean-field result for the n-field model in this regime, xi'f') ~ Gxp(— r/^)/r. 

Finally, for the QC regime, elementary considerations show that the most singular be- 
havior in n(/c) comes from the u = term in the summation over frequency, and hence n(/c) 
is still given by (|2.13| ). This result is intuitively obvious as the QC regime is the interpolation 
regime between the classical and the quantum-disordered ones, and the singularity in n(/c) 
is much stronger in the classical regime. 

C. Gauge field propagator 

In the next sections, we will show that the mean-field (M = oo) behavior of suscepti- 
bilities changes drastically at finite M. The 1/M corrections to susceptibility include the 
propagators of the gauge field. A, and the constraint field, A. We will see that to study 
confinement, it is sufficient to know the propagators of A and A at distances much larger 
then the spin correlation length. At A; ^ the fluctuations of A and of the temporal 
component of the gauge field are well screened P,|lll. However, there is no screening for 



the spatial part of the gauge field p^JTO|JTl[] . The spatial component of the gauge field 
propagator is given by 



9 



/dk — * 
j^Goik, uj)G^{k + q,uj + n)i2k + q)^{2k + g),, (2.15) 

where /cq = = ^- For the RC and the QC regimes, < ksT, and at A; <^ we 
can also restrict to the i7 = component of 11. The evaluation of integrals in the limit of 
q <^ mo is then straightforward, and we find 

n„„„.o)^(*,„-M.)^rE-^. (2.16) 

Inverting now 11^,^ in the Coulomb gauge, we find for the gauge field propagator D^^ at 
small momenta 

DAO) = (^.. - '-f) DiM) = [k. - ^) (2.17) 

where /i(i/) = x, and the values of A are A = ml/kBT (RC regime), A = 2QkBT / 
(QC regime). In the QD regime at T = 0, f2 is a continuous variable, and q,Q are the 
components of a 3D — momentum. We found that eq. (p. 17 ) is still valid in this regime, 
A = 2mo, and fi and u are running over x,y, and r. 

We now consider separately 1/M corrections in the three scaling regimes. 



III. RENORMALIZED-CLASSICAL REGIME 

A. First order in 1/M 

As we said above, we only have to consider the 1/M corrections associated with the 
fluctuations of the gauge field. Corrections related to the fluctuations of the constraint 
field renormalize the spinon mass, but these fluctuations are screened at g < 2mo and are 
therefore irrelevant for the confinement. The first-order 1/M corrections to the polarization 
operator 11 (fc) are shown in Fig. |I|, the propagators and vertex function which appear in 
the 1/M expansion are collected in Fig ^. A simple power counting argument shows that 
the most singular corrections appear if we select a contribution proportional to the external 
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momentum, k, at each interaction vertex with the gauge field. We then obtain near k = 
±2imo, and setting G{p) = G{p,uj = 0), 

ni/M(A:) = - {ksTf J G{p- k/2)G{p + k/2)G{p + I + k/2) 

X {2G{p + k/2) - Gip + I- k/2)) kM^^u - ^f)D{l), (3.1) 

or, in exphcit form. 



x(2 



(p + /) + 5 — i2mo{p cos + Z cos ^z^) p + 5 — i2mQp cos 

X 



+ + 52)2 + 4m2(pcOS(^ + /C0S^)2 (p^ + ^2)2 + 4mg COs2 

^ ' (3.2) 



{{p + /)2 + 52)2 _)_ 4mQ(pcosv9 + Z cosip)'^' 
where ip and t/) are the angles between p and k, and / and k, respectively. The key point in 
the computation of Ili/M{k), as well as in other computations later in the paper, is that for 
6"^ <^ rriQ (which is the only one we consider), the angular integration is confined to a region 
where Lp^ip ±7r/2 |[T^. For practical purposes, it is convenient to restrict the angular 
integration to the vicinity of 7r/2, but extend the integration over p and / from — oo to +cxd. 

Rescaling the angular variables, — > 2^^, ~^ 2m^v then obtain 

' \ 2 

kj^T \ 1 



(3.3) 



"^/"^^^ = -M 1,4^; J^JPJ-J^J-J^J-J'^ip^ + 6^)^ + ^^ 
{p - l)^ + (5^ - i{(p + %l)) p^ + 5"^ - i^p 1 \ 



V ((p-/)2 + 52)2+((^ + V')2 (p2 + 52)2 + <^2 ((^ _ /)2 + ^2)2 + (^^ + ^)2y 

The angular integration extends up to Iv^l, IV'I ~ 0(1/5). The integrals are convergent in 
the ultraviolet, so we can safely set the limits of the angular integration equal to infinity. 
The integration is now trivial, and we find 

We emphasize at this point that only the momenta perpendicular to the external momentum 
k contribute to the angular integration. This implies that the evaluation of the 1/M correc- 
tion to the polarization operator in 2 + dimensions in fact reduces to a one-dimensional 
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problem. This dimensional transmutation will play a key role in our analysis of the se- 
ries of 1/M terms. Finally, substituting D[l) from ( p.l7| ) and performing the momentum 
integration, we obtain 

Notice the absence of divergencies in the integration over / - this is the result of including 
both self-energy and vertex corrections into ni/jv/(^)- 

In obtaining (|3.5| ), we implicitly assumed that 5^ = (/c^ + 4mQ)/4 > 0. Calculations for 
(5^ < proceed along the same lines, and the result is 

\2m?'T 

ni/M(fc)= ^,,,, e + Aml^-^. (3.6) 



M(A;2 + 4m, 



The 1/M corrections to the polarization operator in the RC regime were studied earlier by 
Campostrini and Rossi |jTO|. They numerically evaluated the leading singularity in ni/jv(fc) 



for (5^ +0 and found the same result as in (|3.5|) , with the numerical prefactor 9.425 which, 
as we show, is in fact exactly 37r. However, their estimate for IiiiM{k) for negative 5"^ is 
inconsistent with ( |3.(j| ). 

Let us now discuss the 1/M results. Our first observation is that the 1/5^ divergence 
of liijMik) is stronger than one could expect assuming that n(A;) preserves its form ( [^.13[ ), 
and the two-spinon mass has a regular 1/M expansion (the latter would correspond to 
^i/uik) ~ 0{5~^)). A crude estimate that provides insight into this singularity may be 
obtained by absorbing the 1/M correction into a mass renormalization and solving a self- 
consistent equation for the mass |iy]. One then obtains 11 ~ (A;^ + m^)~^/^, where m = 



2mo(l + 0(1/M^/^)), which is consistent with Witten's result for the mass renormalization. 
However, this self- consistent procedure is clearly not exact, even for large M, because near 
k"^ = — m^, 5 ~ M~^^^, and the 1/M contribution to H has the same order O^l/M"^^^) as 
the leading term. Moreover, simple estimates show that the higher-order 1/M corrections 
behave as 6~^{M6^)~'' and therefore all have the same order near k"^ = —m?. 

We will demonstrate in the next section that the self-consistent procedure |TD[ yielding 
H ~ (A;^+m^)~^/^ is actually incomplete because the 1/M corrections contribute to both the 
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mass renormalization and to the renormalization of the spinon wave function. To see this, 
we will need to solve the problem exactly by summing up the series of the most divergent 
1/M corrections (but we will still assume that M ^ 1 and neglect terms which have relative 
smallness in 1/M). 



B. Solution of ladder series 

The relevant diagrams at each order in 1/M are presented in Fig ^. These diagrams 
form a ladder series. We will first discuss the solution of the ladder series and then show 
that other diagrams are relatively smaller at large M . 

It is convenient to formally include both self-energy and vertex renormalization terms 
into an effective vertex renormalization, such that 

Ii{k) = kBT I ^T{p,k) G{p + k/2)G{p + k/2). (3.7) 

The integral equation for the vertex can be written down as 

r(p, k) = i + ^j r(p + 1 k) G{p + k/2 + i)G{p- k/2 + 1) 

X {2p + k + l)^ {2p-k + I), D^,{1) 

2k T r rl^I 
+^ J 1^ m k)G{p + k/2)G{p + k/2 + f) 

x{2p + k + r)^ {2p + k + l), D^,{1). (3.8) 

Again, the angular integration over both p and / is confined to a narrow region (~ 5) in which 
both internal momenta are nearly orthogonal to the external momentum k. We assume that 
r(p, k) and T{p + /, k) are nonsigular for these values of the angles. We then can set p, / to 
be orthogonal to A; in F and perform angular intergation in the Green functions. Restricting, 
as before, to the integration only near tt/2, and extending the integration over p and / from 
— oo to +00, we obtain 
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where 



Substituting now T{p, k) = ijp' + 5^)\E'(p), we can rewrite ( |3.10| ) as 



(p^ + 5')^>{p) = 1 + (vl/(p) - vl/(p _ /)) . (3.11) 

Finally, performing a Fourier transformation to real space, 

^(^) = r ^e'"-vl/(p), (3.12) 
we obtain for \l/(x) an one- dimensional Schrodinger equation with a source 



+ 'l/(a;)+5'^(a;) = 5(x). (3.13) 



The potential \^(x) is given by 

= !^ f ^Z^(/) (i - e--) = (3.14) 

Notice that the integral over / is free from divergencies. This is again the result of including 
both self-energy and vertex corrections in the ladder series. 
Going back to (|3.9| ), we see that n(A;) takes a simple form 

U(k) = -^^(x = 0). (3.15) 

We emphasize at this point that \I'(x) is not a wave function of the ID Schrodinger equa- 
tion, but rather a Green's function of the inhomogeneous differential equation ( |3.13|) . This 



equation was solved earlier in the content of the weak localization theory , and we simply 
quote the result: 

1 /67rmn\~"^^^ Ai(s) 
^ ^ 2\MJ Ai'{s) ^ ' 

where s = 5'^ (GTrmg/M)^^/^, Ai{s) is the Airy function, and Ai' is its derivative with respect 
to the argument. For positive s, both Ai and Ai' decay exponentially, but for negative s, 
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they have zeros [|T^]. The singular contributions to = 0) indeed come from the zeros of 



Ai'. Expanding near each of these zeros, we obtain 

Ai'{s) ^ (S - Sn) Ai"{Sn) = -\Sn\{s - Sn)Ai{Sn) (3.17) 

where s„ < is the n-th zero of the derivative of the Airy function. Clearly, the polarization 
operator now has a set of simple poles at 5"^ = s„(67rmQ/M)^/^. 

Finally, we remind that we are actually interested in the long-distance behavior of the 
spin correlation function. At long distances, only the pole with the smallest mass is relevant. 
Near this pole, we can approximate the polarization operator and, hence, the staggered 
susceptibility as 

^(*'°> = f (m) FT^- 

where Z= {3n^y/^/\ Si ~ —1.02 [|I8|, and the renormalized mass is given by 



= 4mQ 1 — Si 



'67r\2/3 

m) 

This last result coincides with the expression for the mass obtained by Witten [P 



Eq. ( |3.18D is the key result of this section. We see that the series of divergent 1/M 
corrections near the branch cut not only produce nonanalytical renormalization of the two- 
spinon mass, but also change the branch cut behavior of the staggered susceptibility to a set 
of poles at discrete values of k. In real space, we then have x{r) ~ M~^/^ exp(— mr)/-^? 
at very large distances which agrees with the result for the n— field sigma-model. Notice, 
however, that the residue of each pole contains a factor M~^/^ and vanishes in the limit of 
M = oo. In this limit, the spacing between neighboring poles vanishes, and one recovers 
the mean-field branch-cut solution for n(A;), and hence x(^) ~ exp(— 2mor)/r. It is not 
difficult to show that at finite but large M, the asymptotic behavior associated with the 
lowest pole exists only at very large distances, mr > M^/^, while at smaller distances the 
behavior of spin correlations remains the same as in the mean-field theory. Indeed, at large 
M, the asymptotic behavior is of no practical relevance, since at mr ~ 

Ml/3 > 1, the spm 
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correlation function is already negligibly small. However, for the physical case of M = 2, 
we expect the behavior associated with the pole to dominate at all scales larger than the 
correlation length. 

Notice, however, that the existence of the nonanalytical corrections to the spin correlation 
length makes it very difficult to obtain the exact expression for ^ = in the CP^^^^ model. 
The exact value of ^ in the 0{N) model is known at arbitrary The M = 2 value 

of the overall factor Z in the staggered susceptibility (|3.18|) {Z = 6.517) is also substantially 
larger than Z ~ 2 expected for the 0(3) model from the 1/A^ expansion [||. 

Before concluding this section, we show that the terms not included in the ladder series 
do not contribute to the renormalization of the polarization operator to leading order in 
1/M. The key point is that if we allow the gauge field propagators to cross each other even 
once, the integral over the internal momentum of the two Green functions located between 
the propagators gives zero because the poles are located in the same half-plane. To illustrate 
this, consider the second, "umbrella" -like, diagram in Fig ^. The analytical expression for 
this diagram is 

_ 16m^{kBTY f f r £p(Pl£q ^ sin^Vsin^7 

iiext - 775^ / / / — Z7Z^ — U[L)U[u 



X 



M2 J J J 647r6 KP^ + + 4mgp2 ^os^ 

1 1 



{{p + g)2 + 52)2 + 4mg(]9cosv9 + gcos7) ((p + fy + 52) _^ 2imo(pcosv2 + /cos^^) 

X ^ ^ ^ ^ • (3.20) 

{{p + / + uY + 5^) + 2imo(pcos V9 + / cosip + u cos 7) 

Here ip and 7 are the angles between the external momentum k and the internal ones 
p, Z, and g*, respectively. As before, the relevant M~'^{k'^ + 5^)"^/^ contribution from this 
diagram is confined to the integration over internal momenta which are nearly orthogonal 
to k. Expanding the angles around 7r/2, rescaling them in the same way as in ( p.3| ), and 
shifting the variable ip ^ ip — 'j , we find that the angular integration in ( |3.20| ) reduces to 

1 1 



00 



(i7 / dip 



(p2 + 52)2 + ^2 _ g)2 + ^2)2 + (<^ + ^) 



2 



f°° 1 1 

^ / "^"^((^-02 + 52)+,^ ((p-/-g)2 + <52)+,(y, + ^)- (3.21) 



CO 
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The integration over ip then gives zero. The same is true for all other "crossing" diagrams, 
as well as for the "rainbow" -like graphs for the self-energy (see Fig. Note that this result 
does not depend on the dimensionality. 

IV. QUANTUM-DISORDERED REGION 

We now perform the same type of analysis for the QD regime. This is the low-temperature 
regime at g > Qc, when long-range order in the ground state is destroyed by quantum 
fluctuations. The ground state is then a total spin singlet, and there is a gap. A, towards 
the lowest excited triplet state with S = 1. The temperature corrections to the staggered 
susceptibility in the QD regime are exponentially small, and we therefore restrict our analysis 
in this Section to T = 0. At zero temperature frequency becomes a continuous variable, 
and perturbative 1/M expansion has to be performed in a three-dimensional spacetime. 
The higher dimensionality has already shown up in the computation of the the polarization 
operator at M = oo, which in the QD regime has only a weak, logarithmic singularity at 

A. 1/M corrections 

The computation of the leading 1/M corrections proceeds exactly in the same way as in 
the RC regime. One can easily verify that the typical internal momenta in the 3D analog 
of ( |3.2|) are of the order of 6 and nearly orthogonal to the external momentum k. This, in 
turn, implies that all internal momenta are confined to the plane perpendicular to k, which 
makes the problem effectively two — dimensional. 

Technically this can be seen as follows. Consider 3D vectors p and /. Let them form 
angles and with k, respectively. Then p I = pl{cos if cos + sin ip sin cos 6), where 6 is 
the azimuthal angle (0 < ^ < 27r). In the spirit of our approximation, we expand and ip 
around tt/2. Then immediately p I = pl{cos6 + 0(6'^)) ~ pi cos6, which is the scalar product 
in 2D. 
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Performing the angular integration in the same way as in the previous section, we obtain 

1 r r <Pp<fl D{1) / 1 1 \ 

A-K^M J J 167r4 p2 + ^2 \p2^S^ (p'+f)2 + 52j- ^'^ 

Evaluating this integral, we find 

"^^^ = -4^ (loglA + C), (4.2) 

where C is a constant. As in the RC regime, the 1/M correction to 11 has a stronger depen- 
dence on 5 than is required for a regular 1/M expansion for the two-spinon mass in ( p. 141 ) . 
However, contrary to the previous case, the (log 5) / 5"^ corrections come only from the self- 
energy term. Vertex corrections only contribute to a constant C term in (|4.2|) . Meanwhile, 
if we look back on how the transformation of a branch cut into a pole was obtained in the 
previous section, we can see that this transformation is primarily due to vertex corrections. 
The inclusion of the self-energy terms only serves to make the Fourier transform of D{1) 
infrared- finite, and allows the correct mass renormalization. By analogy, we can expect that 
the series of logarithmic terms will only contribute to the mass renormalization, but the 
mean-field form of the polarization operator will survive. We performed an explicit com- 
putation of the ladder series of the logarithmic 1/M corrections, and found that with the 
logarithmic accuracy, the polarization operator is given by 

n (h\= ^ [ ^^'^ (l l^mglogmp/^ / I2mglogmo/(5 y 
^ ~ Snmo J + S^^ M p^ + \ M p^ + 5^ ) 



1 /■ dp-p 1 ml 

inmo J p2 _|_ ^2 _|_ ^'^^0 mo Ibvrmo k + 



where 



M 



6 , rrin \ . o /, 6 



m' = 4ml[l + — log , ° , = 4m^ ( 1 + -- log M + 0(M) ) . (4.4) 
y M — niQ J \ M / 

We see that niog(A;) still has a logarithmic singularity at = —fn^. However, this result is 
still incomplete, and, in fact, we have to go beyond the logarithmic accuracy and include ver- 
tex corrections, which are crucial for the renormalization of the functional dependence of the 
propagator. Therefore, we must again consider the full ladder equation for the polarization 
operator. 



B. The ladder diagrams 



As before, the ladder equation for the polarization operator is obtained by multiplying the 
equation for the vertex function T{p) (which is a 3D analog of ( p.8|) ) by G{p+ k/2)G{p— k/2) 
and integrating over 3D momentum p. This procedure is completely equivalent to that in 
the RC regime, and we present only the result. Introducing "^{p) = T{p) / {p'^ + 6^), we obtain 

* + = 1 + ^ / ^Dm(p+ T)^'^J ^mmp). (4.5) 

In real space, this equation is again equivalent to the Schrodinger equation with a source 

{~V^+V{r)+S^)^{r)=S{r), (4.6) 
where V{r) is the Fourier transform of the gauge field propagator 

— * 

The integral over / is again infrared finite. The dependence on M in this equation can be 
eliminated by a rescaling r (n^P') ^ ^' obtain 

(^-V% log |f| ^(f) = 5(f), (4.8) 

where the derivatives are taken with respect to f, and we introduced 

e= o U 7-r ^og —-^ ] = 7T(k^ + m^), (4.9) 

where m is given by ([4.4|) . Polarization operator is again related to \E'(r = 0) via n(/c) = 
(7r/2mo)^(r = 0). 

Surprisingly, to the best of our knowledge, no exact solution of the 2D Schrodinger 
equation with a logarithmic potential ("2D Hydrogen atom") is known. However, as 
V(r) ~ log|r| is unbounded from above, there exists a discrete set of energy levels for 
the homogeneous Schrodinger equation, and, just as in the RC regime, ^E'(O) will have an 
infinite number of simple poles at some discrete e„ ~ 0(1) : 
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v[/(0) oc V = ^ V ^" ^ , (4.10) 

Here An is the (unknown) residue of the n-th pole. Using the WKB approximation it is easy 
to estimate that at large n, e„ ~ log(?T,), and An ~ n~^. Also observe that momentum k 
in ( [4.1U| ) is a vector in 2 + 1 space-time dimensions. Hence, eq. ( [4.1U| ) in fact describes the 
dynamic staggered susceptibility. 

Eq. ( 4.101 ) is the key result of this section. We have found that, just as was the case in the 



RC region, the mean-field expression for the staggered susceptibility is also invalid at finite M 
in the QD region, although the corrections to the two-spinon mass and to the susceptibility 
have different dependence on the expansion parameter 1/M in the QD case. The staggered 
susceptibility computed to order 1/M contains an infinite number of simple poles. The 
lowest pole governs the behavior of the spin correlation function at very large distances. 
Performing the Fourier transform of ( 4.10| ), we obtain for equal-time spin-spin correlator 



x{r) ~ (1/M) exp(— mr)/r, compared to the mean-field expression x(r) ~ exp(— 2mor)/r^. 
As in the RC regime, there is a crossover between the two regimes, which occurs at mr ~ M. 

A confinement of spinous due to the logarithmic potential was qualitatively discussed by 
Wen A complimentary scenario of spinon confinement in the QD phase was considered 
by Polyakov [^, and more recently by Read and Sachdev and Murthy and Sachdev PT| . 
They argued that the instanton tunneling events (which we do not consider) lead to a linear 
confining potential between spinous at distances ~ pi ~ 0.06 At large M, 

this scale is much larger than the typical confinement scale, ~ M^, which appears in our 
consideration. The two scales, however, may become comparable at small M. 

V. QUANTUM-CRITICAL REGIME 

Finally, we consider the QC regime, which is the intermediate regime between the RC 
and the QD regimes. At T = 0, the CP^^~^ model action describes D=3 critical theory 
which possess no confinement [^]. The M = oo susceptibility behaves as l/fc, and the 
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1/M corrections from both constraint and gauge fluctuations are of the form \og{l/k)/k. 
Evaluating the 1/M terms and exponentiating the result we obtain x(^) ~ k^^'^~'^'\ where 
Tj = 1 — Notice that the correction to the mean-field result for rj is large. On the 

contrary, for 0{N) sigma- model at the critical point we have rj = which is much closer to 
the Monte-Carlo result for the 0(3) model - r] = 0.03 [Q. 

At finite temperature, the confinement does exist, and we proceed in exactly the same 
way as in the two previous sections. In the QC regime, the typical frequencies are of the 
order of ksT, so, in principle, one has to perform full frequency summation in the 1/M 
formulas for the polarization operator. However, the situation is simpler, as our earlier 
results show that the confining potential is much stronger in the RC regime than in the QD 
regime. As a result, the most divergent contributions in the 1/M series always come from 
the u = terms in the frequency summations. This observation makes the analysis of the 
confinement in the QC regime very similar to that in the RC regime; the only difference is 
in the form of the gauge field propagator. Using ( |2.17| ) and the results of Sec. fT|, we obtain 
after simple manipulations that the branch-cut behavior of the spin susceptibility is indeed 
replaced by a set of poles. Near the lowest pole, the polarization operator behaves as 

where T = (2B^)/v/5 = 0.828471, and the bound state mass is given by 

- = 4eMl-si . (5.2) 



VI. CONCLUSIONS 

We summarize the key results of the paper. We considered here the form of the static 
staggered susceptibility in the CP^^^^ model of an M— component complex unit vector in 
two spatial dimensions. The M = 2 limit of this model describes two-dimensional Heisenberg 
antiferromagnet, and the 1/M expansion we discuss in this paper is a systematic way to 
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calculate various observables. The elementary excitations in the CP^~^ model are S — 1/2 
bosonic spinous. The spin variables are bilinear in spinon fields, and all physical excitations 
are collective modes of two spinous. In the mean-field approximations, the spinous behave 
as noninteracting particles, and the static staggered spin susceptibility has a branch cut 
along the imaginary k axis, terminated at k"^ — — 4mQ, where mo is the mass of a spinon. 

We have shown that in all three scaling regimes at finite T, the fiuctuations of the 
gauge field give rise to divergent 1/M corrections near the branch cut. We selected the 
most divergent corrections and have shown explicitly that they form ladder series. We then 
found that the ladder problem is equivalent to a Schrodinger equation with the 5— functional 
source in dimension D — 1, where D is the spacetime dimension for the original problem. 
Schrodinger equation has a discrete set of solutions, and expanding near each of the solutions, 
we obtained, instead of a branch cut, a sequence of simple poles. This in turn quantitatively 
changes the behavior of the spin correlation function at large distances compared to the 
mean-field predictions. 

The result that the staggered susceptibility has a pole, and not a branch cut, is consistent 
with the results of the alternative perturbative technique for 2D antiferromagnets, which is 
the 1/N expansion for the 0{N) sigma-model for an A?"— component real unit field. The 
latter model describes Heisenberg antiferromagnet at = 3. Elementary excitations in the 
0{N) sigma-model carry S — 1, and, at any N, the static staggered spin susceptibihty has 
a well defined single pole at imaginary momentum k — ±i^~"^, where ^ is the correlation 
length in the system. 

Comparing now for the CP^~^ and the 0{N) sigma-models, we conclude that although 
the differences between them found at the mean-field level are now gone, in practical calcu- 
lations it is more convenient to use the 0{N) model. The reason is that the perturbative 
series in 1/N is regular for this model, and converges much better than the 1/M series for the 
(jpM-i jj^Q^gi ^jig latter can even be nonanalytical in some cases, e.g. when calculating 
the spin correlation length. 

Our results still leave a minor discrepancy between the two approaches: in the 0{N) 

22 



model, the static mean-field susceptibility has a single pole, and corrections do not lead 
to the appearance of new singularities. On the contrary, in the CP^^^^ model, we found, 
to leading order in 1/M, an infinite set of poles. Indeed, 0{N) and CP^^~^ models are 
equivalent only for = 3 and M = 2, when they both describe Heisenberg antiferromagnet. 
However, Bethe-ansatz solution for the 0(3) model also indicates that there is a single mass 
in the problem because the free energy as a function of the uniform magnetic field, /i, has a 
single threshold at he = m [H,^. We therefore expect that all the poles, that we have found 
in the CP^^~^ model, except for the lowest one, should disappear for M = 2. At present, we 
can only speculate how this may occur - the most likely possibility, in our opinion, is that 
the solutions of the Schrodinger equation with n > simply acquire a finite lifetime due to 
higher-order corrections in 1/M. If this is true, then the susceptibility in the CP^^^^ model 
has a single stable pole at any M, while all other poles that we have found to first order in 
1/M, form an incoherent continuum where both real and imaginary part of susceptibility 
are present. This is corroborated by an observation that the n-field propagator of the 0(3) 
model, evaluated at zero frequency and along imaginary k axis, also has an imaginary part 
for \k\ > 3m {k = ik) due to decay of a S* = 1 quanta into three others (the simplest way 
to see this is to substitute a fixed-length constraint on n field by the u{n^Y interaction 



term and do a perturbative expansion in u pj]). We, however, emphasize that this issue is 



irrelevant for the long-distance behavior of the spin correlation function as this behavior is 
associated with the lowest-energy (n = 0) pole, which is always stable. 

A final remark. In this paper, we considered a field-theoretic description of the 
commensurate antiferromagnet. One can also study the anisotropic version of the CP'^'^^^ 
model, with the extra factor (1 — 7) in front of {zd^zl"^ in ( |2.2| ) (0 < 7 < 1). For 7 7^ the 
gauge-field propagator acquires a mass, which prevents the effective potential V{r) in the 
Schrodinger equation from becoming arbitrary large. As a result, even without damping, 
the number of poles in the staggered susceptibility decreases and finally, when 7 exceeds 
some critical value, the disordered phase possesses deconfined S* = 1/2 bosonic spinous []T3[. 



The disordered phase with deconfined spinous was also obtained in by a direct 1/M 
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expansion at small 1 — 7. This limit is of particular interest as the z— field model with 
7 ~ 1 describes incommensurate quantum antiferromagnet near the critical point. 
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FIGURES 



FIG. 1. The first-order 1/M corrections to the polarization operator of spinons. 

FIG. 2. The diagrammatic representation of the spinon and the gauge field propagators and of 
the spinon-gauge field interaction vertex. 

FIG. 3. The full polarization operator and the equation for the full vertex function. The most 
divergent corrections to the mean-field value of Ii{k) near k = ib2imo comes form ladder series. 
For convenience, both vertex and self-energy terms are absorbed into the vertex renormalization. 

FIG. 4. Examples of the second-order diagrams which do not contribute to the ladder series. 
All these diagrams are less divergent near the branch cut singularity than those included in Fig. ^ 
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